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Abstract 

We propose a new decentralized control scheme for DC Islanded microGrids (ImGs) com¬ 
posed by several Distributed Generation Units (DGUs) with a general interconnection topol¬ 
ogy. Each local controller regulates to a reference value the voltage of the Point of Common 
Coupling (PCC) of the corresponding DGU. Notably, off-line control design is conducted in 
a Plug-and-Play (PnP) fashion meaning that (i) the possibility of adding/removing a DGU 
without spoiling stability of the overall ImG is checked through an optimization problem; (ii) 
when a DGU is plugged in or out at most neighbouring DGUs have to update their controllers 
and (iii) the synthesis of a local controller uses only information on the corresponding DGU 
and lines connected to it. This guarantee total scalability of control synthesis as the ImG 
size grows or DGU gets replaced. Yes, under mild approximations of line dynamics, we for¬ 
mally guarantee stability of the overall closed-loop ImG. The performance of the proposed 
controllers is analyzed simulating different scenarios in PSCAD. 
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1 Introduction 


In the recent years, the increasing penetration of renewable energy sources has motivated a growing 
interest for microgrids, energy networks composed by the interconnection of DGUs and loads [1]. 
Microgrids are self-sustained electric systems that can supply local loads even in islanded mode, 
i.e. disconnected from the main grid [2], Besides their use for electrifying remote areas, islands, 
or large buildings, microgrids can be used for improving resilience to faults and power quality in 
power networks [3]. So far, research mainly focused on AC microgrids [1, 2, 3, 4, 5]. However, 
technological advances in power electronics converters have considerably facilitated the operation 
of DC power systems. This, together with the increasing use of DC renewables (e.g. PV panels), 
batteries and loads (e.g. electronic appliances, LEDs and electric vehicles), has triggered a major 
interest in DC microgrids [6, 7, 8]. DC microgrids have also several advantages over their AC 
counterparts. For instance, control of reactive power or unbalanced electric signals are not an 
issue. On the other hand, protection of DC systems is still a challenging problem [8]. 

For AC IrnGs a key issue is to guarantee voltage and frequency stability by controlling inverters 
interfacing energy sources with lines and loads. This problem has received great attention and 
several decentralized control schemes have been proposed, ranging from classic droop control [2, 9], 
to decentralized control [10, 4, 5]. Some control design approaches are scalable, meaning that the 
design of a local controller for a DGU is not based on the knowledge of the whole ImG and the 
complexity of local control design is independent of the ImG size. In addition, the method proposed 
in [4, 5] allows for the seamless plugging-in, unplugging and replacement of DGUs without spoiling 
ImG stability. Control design procedure with these features have been termed PnP [11, 12, 13, 14]. 

Voltage stability is critical also in DC microgrids as they cannot be directly coupled to an 
“infinite-power” source, such as the AC main grid, and therefore they always operate in islanded 
mode. Existing controllers for the stabilization of DC ImGs are mainly based on droop control 
[7, 15]. So far, however, stability of the closed-loop systems has been analyzed only for specific 
ImGs [7, 15]. 

In this paper we develop a totally scalable method for the synthesis of decentralized controllers 
for DC ImGs. We propose a PnP design procedure where the synthesis of a local controller requires 
only the model of the corresponding DGU and the parameters of transmission lines connected to 
it. Importantly, no specific information about any other DGU is needed. Moreover, when a DGU 
is plugged in or out, only DGUs physically connected to it have to retune their local controllers. 
As in [4], we exploit Quasi-Stationary Line (QSL) approximations of line dynamics [16] and use 
structured Lyapunov functions for mapping control design into a Linear Matrix Inequality (LMI) 
problem. This also allows to automatically deny plugging-in/out requests if these operations spoil 
the stability of the ImG. 

In order to validate our results, we run several simulations in PSCAD using realistic models 
of Buck converters and associated filters. As a first test, we consider two radially connected 
DGUs [17] and we show that, in spite of QSL approximations, PnP controllers lead to very good 
performances in terms of voltage tracking and robustness to unknown load dynamics. We also 
show how to embed PnP controllers in a bumpless transfer scheme [18] so as to avoid abrupt 
changes of the control variables due to controller switching. Then, we consider an ImG with 5 
DGUs arranged in a meshed topology including loops and discuss the real-time plugging-in and 
out of a DGU. 

The paper is organized as follows. In Section 2 we present dynamical models of ImGs and 
the adopted line approximation. In Section 3, the procedure for performing PnP operations is 
described. In Section 4 we assess performance of PnP controllers through simulation case studies. 
Section 5 is devoted to some conclusions. 


2 Model of a DC Microgrid 

This section discusses dynamic models of ImGs. For clarity, we start by introducing an ImG 
consisting of two parallel DGUs, then we generalize the model to ImGs composed of N DGUs. 
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Consider the scheme depicted in Figure 1 comprising two DGUs denoted with i and j and connected 
through a DC line with an impedance specified by parameters Rij > 0 and L, 3 > 0. At each DGU 
level, a DC voltage source represents a generic renewable resource and a Buck converter is present 
in order to supply a local DC load connected to the PCC through a series LC filter. For instance, 
the DC load can be a combination of resistive electronic loads and negative resistance of constant 
power loads. Furthermore, we assume that loads are unknown and we treat them as current 
disturbances (Iff) [4, 19]. 


DGU i Line ij and ji DGU j 



Figure 1: Electrical scheme of a DC ImG composed of two radially connected DGUs with unmod¬ 
eled loads. 


Applying Kirchoff’s voltage law and Kirchoff’s current law to the electrical scheme of Figure 
1, it is possible to write the following set of equations: 
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As in [4], we notice that from (lc) and (Id) one gets two opposite line currents I t3 and I 3l . 
This is equivalent to have a reference current entering in each DGU. We exploit the following 
assumption to ensure that Iijit) = — Iji(t ), Vt > 0. 


Assumption 1. Initial states for the line currents fulfill Ii 3 ( 0) = — lji( 0 ). Furthermore, we set 
Li 3 Ij 3 i and dt, 3 — R 3 i . 

Remark 1. According to the terminology in Section 3.f of [20], the system in (lc), (Id) represents 
an expansion of the line model one obtains introducing only a single state variable. System (1) 
can also be viewed as a system of differential-algebraic equations, given by (la)-(lc), (le), (If) 
and lij(t) = -lji(t ). 


At this point, we notice that adopting the above notation for the lines, both DGU models have 
the same structure. In particular, by recalling that the load current II*, * £ i,j is treated as a 
disturbance, (1) is the following linear system 


x(t) = Ax(t) + Bu(t) + Md(t) 
y(t ) = Cx(t) 
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where x = [V), I ti , hj, Iji, Vj t , I t j] T is the state, u = [V ti , V t j] T the input, d = [ILi, lLj] T the 
disturbance and y = [Vi,Vj] T the output of the system. All matrices in (2), which are obtained 
from (1), are given in Appendix A.l. 

Next, we show how to describe each DGU as a dynamical system affected directly by state of 
the other DGU connected to it. An approximate model will be proposed so that there will be no 
need of using the line current in the DGU state equations. 


2.1 QSL model 

As in [16] ai 
QSL model 


As in [16] and [21], we set = 0 and = 0. Consequently, from (lc) and (Id), one gets the 
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By replacing variables lij and Iji in (la) and (le) with the right-hand side of (3), we obtain the 
following model of DGU i 
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Switching indexes i and j in (4) one obtains the model of DGU i 


( x {i] (t) = A u x[q (t) + BiU jij (t) + (t) + (t) 

E w GC/: \ y® (*) = CiX \i\ (*) ( 5 ) 

[ Z[i] (t) = H iy[i] ( t ) 

where x^ = [U:,A,;] T is the state, = V t i the control input, = In the exogenous input and 
z\i] = Vi the controlled variable of the system. Moreover, yu](t) is the measurable output and we 
assume = x\^, while Uu(t) = AijXfjj represents the coupling with DGU j. 

The matrices of r are obtained from (4) and they are provided in Appendix A.2. As 
regards the line, we obtain the subsystem 


E [b] : {"UMllW — Au t ijX[ijj](t) + AnjjX^it) + AijjjXtf](t) 


( 6 ) 


with xu,ij] = lij as the state of the line. The matrices of (6) are derived from (lc) and reported 
in Appendix A.2. We have now all the ingredients to write the model of the overall microgrid 
depicted in Figure 1. In particular, from equations (5) and (6), we get 
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Remark 2. Consider the structure of matrix A 
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We notice that A is block-triangular, therefore its eigenvalues are given by the union of those 

, Aujj and Auji. Moreover, we have Au^j = Auji. By virtue of the positivity 

of the line parameters, line dynamics is asymptotically stable. As a consequence, stability of 
(7) the depends on the stability of local DGUs connected through the QSL model (3). Hence, 
designing decentralized controllers U[*j = At* (j/[*j ) 7 * £ {i,j}, such that the connection of the DGUs 
is asymptotically stable implies stability of the overall closed-loop model of the microgrid. We refer 
to the resulting system as QSL-ImG model. 


of 


4 .. 4 . . 

A A • • 
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2.2 QSL model of a microgrid composed of N DGUs 

In this section, a generalization of model (5) to ImGs composed of N DGUs is presented. Let 
V = (1,... ,A}. First, we call two DGUs neighbours if there is a transmission line connecting 
them. Then, we denote with Mi C D the subset of neighbours of DGU i. We highlight that the 
neighbouring relation is symmetric, consequently j £ Mi implies i £ Mj. In order to describe the 
dynamics of DGU i, we use model (5), with £[,:] = J2je A/i The new matrices of 

are given in Appendix A.3 while the overall QSL-ImG model can be written as follows 


x(t) = Ax(f) + Bu(f) + Md(t) 

X .p'] (t) — An^ijX^pij^ft) A A[jijX pj (t) A A[j : ijX^j^ (f), £ D , Vj £ Mi 


(8a) 

(8b) 


y(f) = Cx(t) 
z (t) = Hy (t) 


(9) 


where x = (xpj,..., X[ N] ) £ R 2N , u = (um,.....,u [JV ]) £ R^, d = (djr],..., d [N] ) £ l w , y = 
(j/[i]> • ■ • 5 V[N]) e R 2N , z = (z m ,... ,Z[ N ]) £ R . Matrices A, Au tij , A Htij , A ljtij , B, M, C and H 
are reported in Appendix A.2 and A.3. 

Note that neither y nor z depend upon states xo^ju Moreover, even x^-] does not influence 
x. Hence, equations (8b) will be omitted in the sequel. 


3 Plug-and-Play decentralized voltage control 

3.1 Decentralized control scheme with integrators 

Let z re f(£) denote the constant desired reference trajectory for the output z (£). In order to 
track asymptotically z re f(£) when d(t) is constant, we consider the augmented ImG model with 
integrators [22]. A necessary condition for having that the steady-state error eft) = z re f (t) — z (t) 
tends to zero as t —> oo, is that for arbitrary constant signals d(f) = d and z re f(£) = z re f, there 
are equilibrium states and inputs x and u verifying 

0 = Ax A Bu A Md 

( 10 ) 
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Proposition 1. 


Given z re f and d, vectors x and u satisfying (11) always exist. 


( 11 ) 
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Proof. From [22], we know that exists x, u verifying (11) if and only if the following two conditions 
are fulfilled: 

(i) The number of controlled variables is not greater than the number of control inputs. 

(ii) The system under control has no invariant zeros (i.e. rank(T) = 3 N). 

Condition (i) is fulfilled since from (5) one has that km and zm have the same size, Vi G T>. In 
order to prove Condition (ii), we exploit the definition of matrices A, B, C and H and the fact 
that electrical parameters are positive. □ 



*[i] 


Z [N] 


Figure 2: Control scheme with integrators for the overall augmented model. 


The dynamics of the integrators is (see Figure 2) 

*[i] W = e [i] W = z ref {i] (t) - Z[H (t) ^ 

— (^) HiCiXfy] (t), 

and hence, the DGU model augmented with integrators is 

! %] (t) = A H X[i\ ( t ) + ( t ) + Mid[i] ( t ) + (t) 

%](*) = C^t) ( 13 ) 

Z[i\ (*) = HiV[i] (*) 


where X[q = [x T n],Vi t } T G M 3 is the state, = i[j] G R 3 is the measurable output, djj] = 
[d[j], z re f^] T G M 2 collects the exogenous signals (both current of the load and reference signals) 

and £[i](t) = AijX[j](t). Matrices in (13) are defined as follows 


An — 


'o o 1 
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Aij — 

i_i 

O' 

0 

Bi = 
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Ci = 

'Ci 
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o' 
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0 1 


Hi = [Hi 0] . (14) 


Through the following proposition we make sure that the pair (An, Bi ) is controllable, thus system 
(13) can be stabilized. 

Proposition 2. The pair ( Au,Bi ) is controllable. 
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Proof. Using the definition of controllability matrix, we get 

M? = [Bi AaBi A\Bi\ 


An 

Bf 


-HiCi 

0_ 



0 Bi AaBi AiBi 
I 0 0 0 


(15) 


Mp. 


Mp „ 


Matrices and Mf' 2 have always full rank, since all electrical parameters are positive, hence 
rank(M G ) = 3. Therefore the pair ( Au,Bi ) is controllable. □ 


The overall augmented system is obtained from (13) as 


{ x(t) = A x(£) + Bu(l) + Md(t) 
y(<) = Cx(t) (16) 

z{t) = Hy(t) 

where x, y and d collect variables and d^ respectively, and matrices A, B, C, M and H 

are obtained from systems (13). 


3.2 Decentralized PnP control 

This section presents the adopted control approach that allows us to design local controlles while 
guaranteeing asymptotic stability for the augmented system (16). Local controllers are synthesized 
in a decentralized fashion permitting PnP operations. Let us equip each DGU Ej?j GC/ with the 
following state-feedback controller 

c [i\ '■ «[»](*) = K iV[i]{t) = K iX[i[it) (17) 

where A',; £ K. lx3 and controllers C^, i £ T> are decentralized since the computation of u^(t) 
requires the state of Sj? GC/ only. Let nominal subsystems be given by S^ GC/ without coupling 
terms £[.;](t). We aim to design local controllers such that the nominal closed-loop subsystem 

= (A u + BiKi)x^(t) + Mid[i](t) 

V[i\(t) = CiX^ft) (18) 

Z[{\(t) = Hiy[i\(t) 

is asymptotically stable. From Lyapunov theory, we know that if there exists a symmetric matrix 
Pi £ K 3 x 3 ,P 2 ; > 0 such that 

(An + BiKi) T Pi + Pi(An + BiKi) < 0, (19) 

then the nominal closed-loop subsystem equipped with controller Cm is asymptotically stable. 
Similarly, the closed-loop QSL-ImG be given by (16) and (17) 

! i(i) = (A + BK)x(f) + Md(t) 

y (t) = C x(t) ( 20 ) 

z(t) = Hy (t) 

is asymptotically stable if matrix P = diag(Pi,..., Pn ) satisfies 

(A + BK) t P + P(A + BK) < 0 (21) 

where A, B and K collect matrices A l; j , B, and K,. for all i, j £V. We want to emphasize that, in 
general, (19) does not imply (21), since one can show that decentralized design of local controllers 
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can fail to guarantee voltage stability of the whole ImG, if coupling among DGUs is neglected (see 
Appendix B in [5] for an example in the case of AC InrGs). In order to derive conditions such 
that (19) guarantees (21), we first define Ad = diag(Ajj,... ,Ann) and Ac = A — Ad- Then, 
we exploit the following assumptions to ensure asymptotic stability of the closed-loop QSL-ImG. 

Assumption 2. (i) Decentralized controllers Cm, i € V are designed such that (19) holds with 


( 22 ) 

where • denotes an arbitrary entry and ??,; > 0 is a local parameter. 

(ii) It holds R r,i Ct . ~ 0, Vi £ V, Vj i £ Mi■ 

As regards Assumption 2-(i), we will show later that checking the existence of P\ as in (22) and 
Ki fulfilling (19) leads to solving a convex optimization problem. On the other hand, there exist 
different ways to fulfill Assumption 2-(ii). In fact, when an upper bound to all ratios R ^ (which 
depend upon line parameters only) is known, one can simply set the control design parameter r]i 
sufficiently small. However, if networks are spread over a small area, the impedances are small 
and predominantly resistive. Therefore, one has R ^ ss 0 by construction and bigger values of 
r]i can be used for synthesizing local controllers. 

Proposition 3. Let Assumption 2 holds. Then, the overall closed-loop QSL-ImG is asymptotically 
stable. 



Proof. We have to show that (21) holds, which is equivalent to prove that 


(A d + BK) P + P(A d + BK) + AgP + PA C < 0. 


(23) 


(“) 


( b ) 


First, we highlight that term (a) is a block diagonal matrix that collects on the diagonal all left 
hand sides of (19). It follows that term (a) is a negative definite matrix. Next, we show that term 
(■ b ) is zero. In particular, each block (i,j) of term ( b ) can be written as 


Pi A-ij 


ApPj 


if j £ Mi 
otherwise 


Using Assumption 2- (ii), we obtain 
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which proves that inequality (23) holds. 


□ 


At this point, in order to complete the design of the local controller Cm, we have to solve the 
following problem. 



Problem 1 . Compute a matrix Ki such that the nominal closed-loop subsystem is asymptotically 
stable and Assumption 2-(i) is verified, i.e. (19) holds for a matrix Pi structured as in ( 22 ). 

Consider the following optimization problem 


O : min an7i + daPi + aaSi 


Yi = 


> 0 


Yi 




< 0 


-Pii gJ 

Gi -I 


< 0 


Yi I 

I Si I 


> 0 

7 i > 0 , j3i > 0 , Si > 0 


( 26 a) 

( 26 b) 

( 26 c) 

( 26 d) 

( 26 e) 


where an, a j 2 and a ,;3 represent positive weights and • are arbitrary entries. Since all constraints 
in (26) are Linear Matrix Inequalities (LMI), the optimization problem is convex and can be solved 
with efficient (i.e. polynomial-time) LMI solvers [23]. 

Lemma 1 . Problem O is feasible if and only if Problem 1 has a solution. Moreover, Ki and P, 
in (19) are given by Ki = GiYf 1 , Pi = Yjf 1 and 11JQ11 2 < y/WA- 

Proof. Inequality (19) is equivalent to the existence of 7 i > 0 such that 

(A u + BJufPi + V, (A u + BiKi) + 7 ” 1 / < 0 (27) 


where Pi is defined in (22). By applying the Schur lemma on (27), we get the following inequality 


(An + BiKi) T Pi + Pi(A u + BiKf) 

I 



(28) 


which is nonlinear in Pi and A,, 
following parametrization trick [23] 


In order to get rid of the nonlinear terms, we perform the 


Yi = p r 1 

Gi = KtYi. 


(29) 


Notice that the structure of Y i: is the same as the structure of Pi. By pre- and post-multiplying 
(28) with and exploiting (29) we obtain 


YiAji + GfBf + AnYi + BiGi Y ' 

Yi - 7 if 


< 0 


(30) 


Constraint (26a) ensures that matrix Pi has the structure required by Assumption 2-(i). At the 
same time, constraint (26b) guarantees stability of the closed-loop subsystem. Further constraints 
appear in Problem O with the aim of bounding ||iCj|| 2 . In particular, we add 1111 2 < s/Jk 
and ||F- || 2 < Si (which via Schur complement, correspond to constraints (26c) and (26d)) to 
prevent ||A',;|| 2 from becoming too large. These bounds imply ||A',;|| 2 < yfflSi and then affect the 
magnitude of control variables. □ 


Next, we discuss the key feature of the proposed decentralized control approach. We first notice 
that constraints in (26) depend upon local fixed matrices ( An , Bi) and local design parameters 
(aii, cti 2 > ai 3 , Pi, Si). It follows that the computation of controller Cm is completely independent 
from the computation of controllers Cu when j 7 ^ i since, provided that problem Vi is feasible, 
controller Cm can be directly obtained through Ki = GfY~ . In addition, it is clear that constraints 
(26c) and (26d) affect only the magnitude of control variables as stated in Lemma 1. Finally, since 
all assumptions in Proposition 3 are verified, the overall closed-loop QSL-ImG is asymptotically 
stable. 
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3.3 Enhancements of local controllers for improving performances 

In the previous section we have shown how to design decentralized controllers Cm guaranteeing 
asymptotic stability for the overall closed-loop system (20). In order to improve transient perfor¬ 
mances of controllers Cm, we enhance them with feed-forward terms for 

(i) pre-filtering reference signals; 

(ii) compensating measurable disturbances. 


3.3.1 Pre-filtering of the reference signal 


Pre-filtering is well known technique used to widen the bandwidth so as to speed up the response 
of the system. Consider the transfer function Fjj](s), from z re f,n(t) to the controlled variable 
Z{i] (t) 


F [t] (s) = (. HA)(sI - {A u + BiKi ))- 1 " 


(31) 


of each nominal closed-loop subsystem (18). By virtue of a feedforward compensator Cm(s), it 
is possible to filter the reference signal z re f ,^ (t) (see Figure 3). Consequently, the new transfer 


Figure 3: Block diagram of closed-loop DGU i with prefilter, 
function from z re f,^(t) to z^(t) becomes 


%s) = q,lWi](s) (32) 

Now, taking a desired transfer function Cm (s) for each subsystem, we can compute, from (32), 
the pre-filter (?[,;] (s) as 

A*]( s ) = -P'[i]( s )- F [*]( s ) _1 (33) 

under the following conditions [22]: 

• Fm (s) must not have Right-Half-Plane (RHP) zeros that would become RHP poles of Cm (s), 
making it unstable; 

• -Fm(s) must not contain a time delay, otherwise C[j](s) would have a predictive action 

• Cm (s) must be realizable, i.e. it must have more poles than zeros. 

Hence, if these conditions are fulfilled, the filter C[j] (s) given by (33) is realizable and asymptoti¬ 
cally stable (this condition is essential since C[i](s) works in open-loop). Furthermore, since F^(s) 
is asymptotically stable (controllers C[j] are, in fact, designed solving the problem Vi), the closed- 
loop system including filters C[.;](s) is asymptotically stable as well. He highlight that, if some 
of the previous conditions are not valid, expression (33) cannot be used. Still, the compensator 
C[j] (s) can be designed for a given bandwidth, as shown in [22]. 
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3.3.2 Compensation of measurable disturbances 


The second enhancement one can introduce regards the compensation of measurable disturbances. 
We remind that, since we assumed that load dynamics is not known, we have modeled the load 
currents for each subsystem of the microgrid as a measurable disturbance dm (t). Let us define 
new local controllers Cm as 

C[{\ : (34) 

Note that Cm are obtained by adding term uw\(t) to the controllers Cm in (17). Hence, (18) can 

! %](t) = (An + BiKi)x[q (t) + Mid[i\(t) + Bi.U{i\(t) 

V\{\ ( t ) = C*%] (t) ■ (35) 

z {i] (t) = Hiy l{ \(t) 


be rewritten as follows 


DGU 


We now use the new 



variable zudt) is 


input U[i\(t) to compensate the measurable disturbance d[,](t) (recall that 
From (35), the transfer function from the disturbance o?m ( t ) to the controlled 


Gi(s) 


(HA)(sI ~ (A u + BiKi ))- 1 


Mi 

0 


(36) 


Moreover, the transfer function from the new input {q.j] (t) to the controlled variable Z[(t) is 


Gi(s) = (HA)(sI - (An + BiKi)ABn (37) 

If we combine (36) and (37), we obtain 

2[i](s) = G i (s)u [i] (s) + Gf(s)d [i] (s). (38) 


In order to zero the effect of the disturbance on the controlled variable, we set 

%]( s ) = N i(s)d[i](s) (39) 

where 

N [{i (s) = -G i (s)- 1 Gf(s) (40) 

is the transfer function of the compensator. Note that 7Vm (s) is well defined under the following 
conditions [22]: 

• G[j](s) must not have RHP zeros that would become RHP poles of 7V[.j](s); 

• G[j](s) must not contain a time delay, otherwise N^(s) would have a predictive action 

• _/V[.j](s) must be realizable, i.e. it must have more poles than zeros. 

In this way, we can ensure that the compensator 7V[j](s) is asymptotically stable, hence preserving 
asymptotic stability of the system. When some of the previous conditions do not hold, formula 
(40) cannot be used and perfect compensation cannot be achieved. Still, the compensator (s) 
can be designed to reject disturbances within a given bandwidth, as shown in [22]. The overall 
control scheme with the addition of the compensators is shown in Figure 4. 

3.4 Algorithm for the design of local controllers 

Algorithm 1 collects the steps of the overall design procedure. 
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Algorithm 1 Design of controller and compensators and A^j for subsystem EjP, GG 
Input: DGU Sg GC/ as in (13) 

Output: Controller Cm and, optionally, pre-filter Cm and compensator A[ji 


(A) Find Kt solving the LMI problem (26). If it is not feasible stop (the controller Cm cannot be 
designed). 

Optional steps 

(B) Design the asymptotically stable local pre-filter Cm and compensator Am as in (40). 


3.5 PnP operations 

In the following section, the operations for updating the controllers when DGUs are added to or 
removed from an ImG are presented. We remind that all these operations are performed with the 
aim of preserving stability of the new closed-loop system. Consider, as a starting point, a microgrid 
composed of subsystems EjP GC/ ,« € V equipped with local controllers C^j and compensators C^j 
and Am, i £V produced by Algorithm 1. 

Remark 3. In order to avoid jumps in the control variable when local regulator are switched, we 
embedded each local regulator into a bumpless control scheme [18] shown in Appendix B. 

Plugging-in operation Assume that the plug-in of a new DGU Ej^_ G j described by ma¬ 
trices, Ajv+ijv+i, Bjv+i, Cjv+i, Mjv+i) Hn+i and (A.v+i j}jeU N +i needs to be performed. Let 
Mn+i be the set of DGUs that are directly coupled to Ef^b through transmission lines and let 
{An + i j}j G jg- N+1 be the matrices containing the corresponding coupling terms. According to our 
method, the design of controller C[jv+i] and compensators Crjv+i] and Apv+i] requires Algorithm 
1 to be executed. Since DGUs E GGC/ , j G A/jv+i, have the new neighbour E^_U, we need to re¬ 
design controllers Cm and compensators Cy i and Am, Vj £ Nn+i because matrices Ajj, j € Wjv+ 1 
change. 
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Only if Algorithm 1 does not stop in Step A when computing controllers Cry for all k G 
A/jv+i U {N + 1}, we have that the plug-in of Eis allowed. Moreover, we stress that the 
redesign is not propagated further in the network and therefore the asymptotic stability of the 
new overall closed-loop QSL-ImG model is preserved even without changing controllers Cm, Cm 
and JV W , * ^ {./V A 1} U x. 

Prior to real-time plugging-in operation (hot plugging-in), it is recommended to keep set points 
constant for a sufficient amount of time so as to guarantee control variable in the bumpless con¬ 
trol scheme (see Remark 3) is in steady state. This ensures smooth behaviours of the electrical 
variables. 

Unplugging operation Let us now examine the unplugging of DGU Ej^ G!7 , k £ T>. The 
disconnection of Ej^ G,7 from the network leads to a change in matrix Ajj of each S GG!7 , j £ A4- 

Consequently, for each j £ A4, we have to redesign controllers Cui and compensators Cm and Afoi. 
As for the plug-in operation, we run Algorithm 1. If all operations can be successfully terminated, 
then the unplugging of E^ GC/ is allowed and stability is preserved without redesigning the local 
controllers Cm, j ^ A4- 

When an unplugging operation is scheduled in advance, it is advisable to follow an hot unplug¬ 
ging protocol similar to the one introduced for the plugging-in operation. 


4 Simulation results 

In this section, we study performance brought about by PnP controllers described in Section 3. 
As a starting point, we consider the ImG depicted in Figure 1 with only two DGUs (Scenario 1) 
and we evaluate performance in terms of tracking step references as well as hot plugging-in of the 
two DGUs and robustness to unknown load dynamic. Then, we extend the analysis to an ImG 
with 6 DGUs (Scenario 2) and we show that stability of the whole microgrid is guaranteed. 

Simulations have been performed in PSCAD, a simulation environment for electric systems 
which allows to implement the microgrid model with realistic electric components. 

For both scenarios, we run a simulation from time 0 s up to time 10 s. Each simulation has 
been split into subparts that are discussed next. 

4.1 Scenario 1 

In this Scenario, we consider the ImG shown in Figure 1 composed of two DC DGUs connected 
through high resistive-inductive lines supporting 10 and 6 O loads, respectively. For the sake of 
simplicity, we set i = 1 and j = 2. The output voltage reference v* MG has been selected at 48 V 
and it is equal for both DGUs. Parameters values for all DGUs are given in Table 1 in Appendix 
C. Notice that that they are comparable to those used in [17] and [7]. 

4.1.1 Voltage reference tracking at the startup 

We assume that at the beginning of the simulation (t = Os), subsystems E?, GC/ and Ep, GC/ are 
not interconnected. Therefore, stabilizing controllers Ci, i = 1,2 are designed neglecting coupling 
among DGUs. Moreover, in order to widen the bandwidth of each closed-loop subsystem, we use 
local pre-filters C[ j], i = l,2of reference signals. The desired closed-loop transfer functions F)(.s), 
* = 1,2 have been chosen as low-pass filters with DC gain equal to 0 dB and bandwidth equal to 
100 Hz. The eigenvalues of the two decoupled closed-loop QSL subsystems are shown in Figure 5a. 
Moreover, by running Step B of Algorithm 1 we obtain two asymptotically stable local pre-filters 
Ci, i = 1,2 whose Bode magnitude plots are depicted in Figure 5b. Notice that through the 
addition of the pre-filters, the frequency response of the two closed-loop transfer functions Fi(s), 
i = 1,2 coincide with the frequency response of the desired transfer functions Fi(s), i = 1,2 (see 
the green line in Figure 5c). Figures 6a and 6b show the voltages at PCC\ and PCC 2 - Note that 
the controllers ensure an excellent tracking of the reference signals at the startup in a very short 
time. 
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Eigenvalues Closed-Loop (CL) 



(a) Eigenvalues of the two decoupled closed-loop QSL 
subsystems. 



(b) Bode magnitude plot of pre-filters Cm, i = 1, 2. 



(c) Bode magnitude plot of Fi(s), i = 1, 2 with (green) 
and without (blue) pre-filters. 

Figure 5: Features of PnP controllers for Scenario 1 when the DGUs are not interconnected. 




(a) Load voltage at PCC \. (b) Load voltage at PC( '2 ■ 

Figure 6: Scenario 1 - Voltage reference tracking at the startup. 


4.1.2 Hot plugging-in of DGUs 1 and 2 

At time t = 2 s, we connect DGUs 1 and 2 together. This requires real-time switch of the local 
controllers which translates into two hot plugging-in operations as described in Section 3.5. The 
new decentralized controllers for subsystems E GGu and E^ GC/ are designed running Algorithm 1. 
Notice that the interconnection of the two subsystems lead to a variation of each DGU dynamics, 
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therefore even compensators Cm and AU, i = 1,2 need to be updated. In particular, the new 
desired closed-loop transfer functions Fi(s), i = 1,2 have been chosen as low-pass filters with DC 
gain equal to 0 dB and bandwidth equal to 100 Hz. 

Since Algorithm 1 never stops in Step A, the hot plug-in of the DGUs is allowed and local 
controllers get replaced by the new ones at t = 2 s. Figure 7a shows the closed-loop eigenvalues 
of the overall QSL ImG composed of two interconnected DGUs. The Bode magnitude plots of 
compensators CCi and AUi, * = 1,2 are depicted in Figure 7b and 7c, respectively, while the 
singular values of the overall closed-loop transfer function F(s) with inputs [z re f m , z re f [2] ] T and 
outputs [ 2 [i], Z[ 2 ]] T are shown in Figure 7d. 


Eigenvalues Closed-Loop (CL) 




(a) Eigenvalues of the closed-loop QSL microgrid with 
(red) and without (blue) couplings. 


(b) Bode magnitude plot of pre-filters Cm, i = 1, 2. 




(c) Bode magnitude plot of disturbances compensators (d) Singular values of F(s) with (green) and without 
, i = 1,2. (blue) pre-filters. 

Figure 7: Features of PnP controllers for Scenario 1 when the DGUs are connected together. 


Figure 8 shows the dynamic responses of the voltages at PCC\ and PCCi when the subsys¬ 
tems are connected together. We highlight that the bumpless control transfer schemes ensure no 
significant deviations in the output signals when the controller switch is performed. Moreover, 
through the proposed decentralized control strategy, voltage regulation is excellent. 


4.1.3 Robustness to unknown load dynamics 

Next, we assess the performance of PnP controllers when loads suddenly change at a certain time. 
To this purpose, at t = 3 s we decrease the load resistances at PCCi and PCC 2 to half of their 
initial values. Figure 9 shows the response of the ImGs. Figures 9a and 9b show the load voltage 
at PCCi and PCC 2 which confirm very good compensation of the current disturbances produced 
by load changes. We notice small oscillations of the voltage signals due to the presence of complex 
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48.2 


48.2 


48.1 

48 

47.9 

47.8 


2.2 2.4 2.6 2.8 

time (s) 

(a) Load voltage at PCC \. 


48.1 

48 

I 

47.9 

47.8 


2.2 2.4 2.6 2.8 

time (s) 

(b) Load voltage at PCC 2 - 


Figure 8: Scenario 1 - Impact of bumpless control transfer on the hot plug-in at time t = 2 s. 


coniugate poles in the transfer function of the closed-loop overall system including couplings (as 
shown in Figure 7a). However, these oscillations disappear after a short transient. We recall 
that load currents (see Figures 9c and 9d) are treated as measurable disturbances in our model. 
Varying the load resistance, induces step-like changes in the disturbances. 



(a) Voltage at PCC \. 




(b) Voltage at PCC 2 . 



(c) Instantaneous load current II i. 


(d) Instantaneous load current Il 2 - 


Figure 9: Scenario 1 - Performance of PnP decentralized voltage control in presence of load switches 
at time t = 3 s. 


4.1.4 Voltage tracking for DGU 1 

Finally, we evaluate the performance in tracking step changes in the voltage reference at one PCC 
(e.g. PCC\) when the DGUs are connected together. This test is of particular concern if we 
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look at the concrete implementation of islanded DC microgrids. In fact, changes in the voltage 
references can be required in order to regulate power flow among the DGUs, or to control the 
state-of-charge of possible batteries embedded in the ImG. 

To this purpose, at t = 4 s we let the reference signal of DGU 1, u* MG , step down to 47.6 V. 
Notice that this small variation of the voltage reference at PCC\ is sufficient to let an appreciable 
amount of current flow through the line, since the line impedance is quite small. The dynamic 
responses of the overall microgrid to this change are shown in Figure 10. As one can see, controllers 
guarantees good tracking performances in a reasonable time with small interactions between the 
two DGUs. 




(a) Voltage at PCC\. (b) Voltage at PCC' 2 ■ 

Figure 10: Scenario 1 - Performance of PnP decentralized voltage controllers in terms of set-point 
tracking for DGU 1. 


4.2 Scenario 2 

In this second scenario, we consider the meshed ImG shown in Figure 11. As one can notice, the 
main difference with respect to Scenario 1 is that some DGUs have more than one neighbour. This 
means that the disturbances influencing their dynamics will be greater. Moreover, the presence of 
a loop further complicates voltage regulation. To our knowledge, control of loop-interconnected 
DGUs has never been investigated for DC microgrids. 



Figure 11: Scenario 2 - Scheme of the ImG composed of 5 DGUs, connected through transmission 
lines (black arrows). 

In order to assess the capability of the proposed decentralized approach to cope with heteroge¬ 
neous dynamics, we consider an ImG composed of DGUs with non-identical electrical parameters. 
They are listed in Tables 2, 3 and 4 in Appendix C. 

We also assume that DGUs 1-5 supply 10 D, 6 11, 20 D, 2D and 4 D loads, respectively. 
Moreover, we highlight that, for this Scenario, no compensators C,; and iVj have been used. At the 
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beginning of the simulation, all the DGUs are assumed to be isolated and not connected to each 
other. However, we choose to equip each subsystem Ti GGU , i £ V = {1,..., 5}, with controllers 
Cy j] designed by running Algorithm 1 and taking into account couplings among DGUs. This is 
possible because, as shown in Section 3.2, local controllers stabilize the ImG also in absence of 
couplings. Because of this choice of local controllers in the startup phase, when the five subsystems 
are connected together at time t = 1.5 s, no bumpless control scheme is required since no real-time 
switch of controllers is performed. The closed-loop eigenvalues of the overall QSL ImG are depicted 
in Figure 12a while Figure 12b shows the closed-loop transfer function of the whole microgrid. 


Eigenvalues Closed-Loop (CL) 


600 

400 

200 

0 * 

-200 

-400 

-600 


x eig CL with coup ings 
. eig CL without couplings 


C OnC 0 





(a) Eigenvalues of the closed-loop QSL microgrid with , Q . . . , , , 

(red) and without (blue) couplings. ' mgu ar va ues o (s) 


Figure 12: Features of PnP controllers for Scenario 2 with 5 interconnected DGUs. 


4.2.1 Plug-in of a new DGU 

For evaluating the PnP capabilities of our control approach, we simulate the connection of DGU 
with Sj°] GC/ and S^ GC/ , as shown in Figure 13. Therefore, we have A/e = {1,5}. In 



Figure 13: Scenario 2 - Scheme of the ImG composed of 6 DGUs after the plugging-in of T,^ GU . 

principle, subsystems , j € Afe must update their controllers Cyjy (see Section 3.5). However, 

we highlight that previous controllers for DGUs ^y^F u and ^y^ GU can be also maintained, provided 
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that the already computed matrices Kj , j £ Afe still fulfill all constraints in (26) for the new ImG 
topology. Since this test succeeds, we proceed by executing Algorithm 1 for synthesizing C[ 6 ] for 
the new DGU only. Algorithm 1 never stops in Step A and therefore the addition of is 

allowed. The real-time plugging-in of Ej^ GC/ is executed at time t = 2 s. Until the plug-in of 
Gc/ , common reference v* MG for DGUs 1-5 is the same as for DGUs 1-2 in Scenario 1 and the 

subsystem E^ GG is assumed to work isolated, tracking the reference voltage v* MG . Figures 14a 
and 14b show respectively the closed-loop eigenvalues and the singular values of the closed-loop 
F{s) of the overall QSL ImG in Figure 13 equipped with the controllers described above. From 
Figure 15, we note that right after the hot plug-in of E^? GC/ at t = 2 s, load voltages of Ej^ GC/ and 

E GGC/ do not deviate from the respective reference signals. 



(a) Eigenvalues of the closed-loop QSL microgrid with . . , 

(red) and without (blue) couplings. ln S u ar va ues o (s) 


Figure 14: Features of PnP controllers for Scenario 2 with 6 interconnected DGUs 


4.2.2 Robustness to unknown load dynamics 

In order to test the robustness of the overall ImG to unknown load dynamics, at t = 3 s we halve 
the load of DGU 6 , which was equal to 8 D for t < 3 s. 

Figures 16a and 16b show that, when the load change of Ej^ GC/ occurs, the voltages at PCC\ 

and PCC5 exhibit very small variations which last for a short time. Then, load voltages of E GGC/ 
and E GGC/ converge to their reference values. Similar remarks can be done for the new DGU 

E^ gc/ : as shown in Figure 16c, there is a short transient at the time of the load change, that is 
effectively compensated by the control action. These experiments highlight that controllers Cm, 
i = 1 ,..., 6 may ensure very good tracking of the reference signal and robustness to unknown load 
dynamics even without using compensators C[6] and _/V[ 6 ]. 

4.2.3 Unplugging of a DGU 

Next, we simulate the disconnection of E^ GC/ so that the considered ImG assumes the topology 
shown in Figure 17. The set of neighbours of DGU 3 is A 3 = {1,4}. 

Because of the disconnection, there is a change in the local dynamics Ajj of DGUs E GGG , 
j € A 3 . Then, in theory, each controller Cyj, j € A 3 must be redesigned (see Section 3.5). As for 
the plugging-in operation, we decide to maintain the previous controller for DGUs 1 and 4, after 
checking that the already computed matrices Kj, j £ A 3 fulfill all constraints in (26) even when 
DGU 3 is removed. Since this test ends successfully, the disconnection of Ep GC/ is allowed. Figure 
18a shows that the closed-loop model of the new QSL microgrid is still asymptotically stable in 
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(b) Voltage at PCC 5 . 
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(c) Voltage at PCCq. 

Figure 15: Scenario 2 - Performance of PnP decentralized voltage controllers during the hot plug-in 
of DGU 6 at time t = 2 s. 


spite of the unplugging operation while Figure 18b shows the closed-loop transfer function F(s) of 
the ImG. Hot-unplugging of Epj 3u is performed at time t = 7 s. As shown in Figure 19, the load 

voltages of DGU E^ GC/ , j £ A /3 do not deviate from the respective reference signals. We stress 
again that stability of the microgrid is preserved despite the disconnection of E^ GC/ . 

5 Conclusions 

In this paper, a decentralized control scheme for guaranteeing voltage stability in DC ImGs was 
presented. The main feature of the proposed approach is that, whenever a plugging-in or -out 
of DGUs is required, only a limited number of local controllers must be updated. Moreover, as 
mentioned in Section 4.1.4, local voltage controllers should be coupled with a higher control layer 
devoted to power flow regulation so as to orchestrate mutual help among DGUs. To this purpose, 
we will study if and how ideas from secondary control of ImGs [7] can be reappraised in our 
context. 
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(b) Voltage at PCC 5 . 



Figure 16: Scenario 2 - Performance of PnP decentralized voltage controllers in terms of robustness 
to an abrupt change of load resistances at time t = 3 s. 



Figure 17: Scenario 2 - Scheme of the ImG composed of 5 DGUs after the unplugging of S^ GC/ . 
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Eigenvalues Closed-Loop (CL) 




(a) Eigenvalues of the closed-loop QSL microgrid with . > . . 

(red) and without (blue) couplings. ' ingu ar va ues o (s). 

Figure 18: Features of PnP controllers for Scenario 2 after the unplugging of DGU 3. 




Figure 19: Scenario 2 - Performance of PnP decentralized voltage controllers during the hot- 
unplugging of DGU 3 at t = 7 s. 
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A Matrices appearing in microgrid models 

The following appendix collects all matrices appearing in Section 2. 


A.l Matrices in the model (2) 


A = 


Li 

1 

Lji 

0 

0 


1 

c » 

_ t\,tz 

Rti 

0 

0 

0 

0 


1 

Cti 

0 

0 

0 

0 

0 

_ Rjj 

0 

1 

Lij 

Lij 

0 

_ Rjj 

_ 1 _ 

Lji 

Lji 

0 

1 

Ctj 

0 

0 

0 

_ 1 _ 

Ltj 


0 

0 

0 

0 

1 

Ctj 
_ An 

Ltj 


B = 


' 0 

0 ■ 


'1 

o' 


r 1 

Cti 

0 ' 

1 

Lti 

0 


0 

0 


0 

0 

0 

0 

0 

0 

C T = 

0 

0 

M = 

0 

0 

0 

0 

0 

0 

o o 

°H4 


o o 

1 

0 


o o 

_1 

Ctj 

0 


A.2 Matrices in the QSL model (5) and (6) 
DGU-z, i e {1,2} 


T i ; — 


Rij Cti 


Rt, 
L u 


Aij — 


Rij Cti 

0 



' o " 


■_ 1_- 


'l 0’ 

Bi = 

1 

. Lti - 

Mi = 

Cti 

0 

Ci = 

0 1 


Hi = [1 0 ] 


Line i ^ j 


A i , j j — 




A _ 

~ r - 

Lji 


(41) 


A.3 QSL model of microgrid composed of N DGUs 


DGU-i, i e V 


A-a — 


)J eA fi Rij Cti 

_ ]_ 

Lti 


1 


C ti 


Ru 

Lti 


Aij — 


1 

Rij Cti 

0 


0 

0 


(42) 


(43) 


We remind that Rij and are the resistance and the inductance of the line between DGU i and 
DGU j. Moreover, matrices Bj, Cj , M t and Hi are equal to those appearing in Section A.2. 
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Overall model of a microgrid composed by N DGUs 
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B Bumpless control transfer 

Since the controller is a dynamic system, it is necessary to make sure that the state of the system 
is correct when a switch of the controller (i.e. a plugging-in or unplugging operation) is required. 
Assuming that the control switch is made at a certain point in time t, we call u pre c,i the control 
signal produced by the controller Ci up to time t. It might happen that the updated controller 
will provide a control variable Ui different from u prec j j. Therefore, it is necessary to ensure there 
is no substantial change in the two outputs at t. This is called bumpless control transfer [18]. 

A bumpless control transfer implementation of PnP local controller for system T, GGU is illus¬ 
trated in Figure 20. 


C (s) 

prec' ' 


C(s) 






r(s) K> 


|—*V« 



r(s)-s 


r(s) 


N(s) 




XT. 


Figure 20: Bumpless control transfer implementation. 


For the sake of simplicity we drop here the index % of the subsystem and associated local 
variables and all switches in Figure 20 are assumed to commute at time t. In the Figure, vector 

K = [k v k c ki] T 


contains the parameters of the local controller to be activated at time t. Notice that the integrator 
embedded in the DGU model for zeroing the steady-state error is replaced by block A (highlighted 
in red in Figure 20), where the polynomial T(s) has to be chosen such that ki > T(0) and the 
transfer function 


F(s) - 




T(s 


is asymptotically stable and realizable. In block A, a switch is present so that the signal is either 
Uprec (up to time t) or u (right after t). The variable u pre c is given by 


Uprec — Uprec k v V k c I± U (45) 

where u is the additional input produced by compensator N(s), computed with respect to the 
dynamics of the system after the commutation (set N(s ) = 0 if such a compensation is not 
implemented). We highlight that since there could be a transient in the u response to track signal 
Uprec , it is fundamental to wait for the two signals to become similar before proceeding with the 
commutation. In this way we avoid jumps in the control variable. 

Furthermore, if an optional prefilter of the reference is implemented, at time t is also necessary 
to commute from transfer function C pre c(s) to C(s), since each plugging-in or unplugging operation 
of other DGUs in the overall ImG lead to a variation of the local dynamics of the considered 
subsystem T, GGU (see the term ~cTi in ( 42 ))- 
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C Electrical and simulation parameters of Scenario 1 and 2 

In this appendix, we provide all the electrical and simulation parameters of Scenarios 1 and 2 
(which are described in Sections 4.1 and 4.2, respectively). 


Parameter 

Symbol 

Value 

DC power supply 

Vdc 

100 V 

Output capacitance 

c u 

2.2 mF 

Converter inductance 

Lt* 

1.8 mH 

Inductor + switch loss resistance 

Rt* 

0.2 n 

Switching frequency 

fsw 

10 kHz 

Transmissione line inductance 

L*o 

1.8 

Transmission line resistance 

R*o 

0.05 n 


Table 1: Electrical setup and line parameters 


Table 2: VSC filter parameters for DGUs T. GGU , i 6} in Scenario 2. 


DGU 

Resistance 

Capacitance C t (m 

F ) Inductance L t {mH) 

S'DGU 

w 

0.2 

2.2 

1.8 

k'DGU 

I 2 ! 

0.3 

1.9 

2.0 

X'DGU 

[3] 

0.1 

1.7 

2.2 

x^DGU 

M 

0.5 

2.5 

3.0 

x^DGU 

^[5] 

0.4 

2.0 

1.2 

X'DGU 

^[6] 

0.6 

3.0 

2.5 

Table 3: Transmission lines parameters for Scenario 2. 

Connected DGUs (i, j) 

Resistance i? s (D) 

Inductance L s (/il7) 


(1,2) 

0.05 

2.1 


(1,3) 

0.07 

1.8 


(3,4) 

0.06 

1.0 


(2,4) 

0.04 

2.3 


(4,5) 

0.08 

1.8 


(1,6) 

0.1 

2.5 


(5,6) 

0.08 

3.0 
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Table 4: Common parameters of DGUs SP, GC/ , i = {1,..., 6} in Scenario 2. 


Parameter 

Symbol 

Value 

Electrical parameters 

DC power supply 

Vdc 

100 V 

Output capacitance 

c t 

2.2 mF 

Converter inductance 

L t 

1.8 mH 

Inductor+switclr loss resistance 

Rt 

0.2 n 

Switching frequency 

fsw 

10 kHz 
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